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(~| ■ Abstract 

Q^' We present the analytical description of top-antitop pair production near the thresh- 

' old in e+e^ annihilation and 77 collisions. A set of basic observables considered in- 

* * ■ 

> ■ 

, eludes the total cross sections, forward-backward asymmetry and top quark polariza- 

X: 

H tion. The threshold effects relevant for the basic observables are described by three 

universal functions related to S wave production, P wave production and S — P inter- 
ference. These functions are computed analytically up to the next-to-next-to-leading 
order of NRQCD. The total e^e~ tt cross section near the threshold is obtained 
in the next-to-next-to-leading order in the closed form including the contribution due 
to the axial coupling of top quark and mediated by the Z-boson. The effects of the 
running of the strong coupling constant and of the finite top quark width are taken 
into account analytically for the P wave production and S — P wave interference. 
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1 Introduction. 



Being heavy the top quark undergoes fast weak decays. The relatively large width Tt of the 
top quark is mainly saturated by the decay channel t — >■ Wb and keeps the effective energy 
of top-antitop system in the complex plane far enough from the cut along the positive 
semiaxis. Thus it serves as a sufficient infrared cutoff for long distance effects avoiding the 
problem of strong coupling. This allows one to bypass possible nonperturbative regions and 
is the key observation for the theoretical study of the top-antitop pair production near the 
two-particle threshold Because the relevant scale i/Ftrnt, with rrit being the top quark 
mass, is much larger than Aqcd, the QCD perturbation theory expansion is applicable for 
the theoretical description of physical phenomena near the top quark threshold if singular 
Coulomb effects are properly taken into account [|l], 0, This feature turns the processes 
involving the top quarks into a unique laboratory for perturbative investigation of threshold 
effects. Experimental study of the top-antitop pair threshold production is planned to be 
performed at the Next Linear Collider both in high energy e~^e~ annihilation and 77 collision 
1^. High quality experimental data that can be obtained in such experiments along with a 
very accurate theoretical description of them make the processes of top-antitop pair threshold 
production a promising place for investigating quark-gluon interactions. This investigation 
concerns both general features of interaction and precise quantitative properties such as the 
determination of numerical values of the strong coupling constant a^, the top quark mass, 
and the top quark width. Though the main features in both e+e~ and 77 processes of top 
quark pair threshold production are rather similar the strong interaction corrections and 
relativistic corrections are different for them. Therefore a simultaneous analysis of these two 
processes extends possibilities of studying fine details of the top quark threshold dynamics. 
Besides the total cross sections which are mainly saturated by the 5* wave final state of the 
top quark-antiquark pair, there is a set of observables sensitive to the P wave component. 
For example, the 5* and P partial waves of the final state top quark-antiquark pair produced 
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in 77 collisions can be separated by choosing the same or opposite helicities of the colliding 
photons 0. This gives an opportunity of direct measurement of the P wave amplitude 
which is strongly suppressed in the threshold region in comparison with the S wave one. On 
the other hand, the forward-backward asymmetry of the quark-antiquark pair production in 
e~^e~ annihilation || and top quark polarization [|, 0] are determined by the S — P wave 
interference in both processes. This provides us with two additional independent probes of 
the top quark interactions. 

The finite order perturbation theory of QCD breaks down in the threshold region of par- 
ticle production due to the presence of singular {as/ (3)"' Coulomb terms. Here /3 is a velocity 
of the heavy quark. However the resummation of these Coulomb contributions which are 
most important quantitatively in the threshold region is possible and can be systematically 
done in the framework of nonrelativistic QCD (NRQCD) (for the recent development 
of the NRQCD effective theory approach see [|, |10|, |Tl|, 0, |13|, 0, |T6l). Note that the 
characteristic scale of the Coulomb effects for the top quark production asfrit is compara- 



ble numerically with the cutoff scale for infrared effects y/Ttrrit so the Coulomb effects are 
not suppressed by the top quark width. The determination of higher order corrections in 
the QCD coupling constant and of relativistic corrections in the case when Coulomb ef- 
fects have to be taken into account beyond the finite order perturbation theory requires 
the perturbative expansion for the complete correlator to be performed near the Coulomb 
approximation rather than near Green functions of the free theory which is the standard 
pattern of perturbation theory for the infrared safe high energy processes. 

Recently has been made a rather essential progress in the theoretical description of heavy 
quark-antiquark threshold dynamics within NRQCD. The evaluation of next-to-leading order 
(NLO) and next-to-next-to-leading order (NNLO) corrections to the heavy quark threshold 
production in e~^e~ annihilation has been done both within the analytical approach fl^, |18 



|T|, 0, 0, |2|, |2|] and numerically |2|, g^, ^, |9|, while the NLO corrections 
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to the heavy quark threshold production in 77 colhsion where computed analytically 
The analysis of NNLO corrections in the last case is still absent. However, this analysis is 
necessary for the accurate quantitative study of the process since the NNLO contribution 
is found to be relatively large in the case of the top quark production in e^e~ annihilation 
p6| , ^ and one can expect that some large corrections emerge also in the case of the top 
quark threshold production in 77 collision. Moreover, a semianalytical analysis of the high 
order corrections to the top quark threshold production cross section in e^e~ annihilation 
has been performed so far I^B], ^ while the essential part of corrections has been 

accounted for numerically . Therefore the complete analytical description of the process is 
also desirable^. Furthermore, the forward-backward asymmetry and top quark polarization 
has been analyzed in NLO only numerically as well as the axial contribution in the 
e"^e^ ti process p8|. In this case the numerical study is more involved because of necessity 
to construct the Green function for the P wave which leads to the more singular differential 
equations in comparison with the S wave. The case of P wave production in 77 collision 



0, ^ clearly demonstrates that the numerical analysis with an explicit cutoff of the 
hard momentum contribution is not sufficient for an accurate account for the finite top quark 
width for these quantities because the relativistic effects are not taken into account properly. 

In the present paper we give a simultaneous analysis of several observables relevant to 
e+e^ ti annihilation and 77 ti collisions near the top quark production threshold in 
high orders of NRQCD. The total cross sections are computed in NNLO of NRQCD which 
includes a^, and corrections in the coupling constant and the heavy quark velocity 
f3 to the nonrelativistic Coulomb approximation. Explicit analytical expressions for the soft 
part of corrections are obtained. The threshold cross section of the ti production in e+e~ 
annihilation is obtained in the closed form including the contribution due to the top quark 
axial coupling. The hard part of the correction to the 77 ti threshold cross section is 



^When this work was in its final stage a letter appeared where the photon mediated top quark 
production in e+e^ annihilation was analyzed analytically. 
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found with the logarithmic accuracy. The forward-backward asymmetry of the top quark- 
antiquark pair production in e'^e~ annihilation and top quark polarization in both processes 
of e+e~ annihilation and 77 collisions are computed up to NLO. 

The paper is organized as follows. In the next Section the nonrelativistic approximation 
for the basic observables of top quark- ant iquark pair production near the threshold is formu- 
lated. In Section 3 the threshold effects are described by three universal functions related to 
the S, P wave production and S — P wave interference which have been computed analyti- 
cally within NRQCD. In Section 4 we present our numerical analysis and the discussion of 
the obtained results. The last Section is devoted to our conclusions. Some explicit analytical 
formulae are given in Appendix. 

2 The nonrelativistic approximation near the produc- 
tion threshold. 

In this section we describe the set of observable which will be analyzed: the total cross 
sections, the forward-backward asymmetry, the polarization of top quark. We formulate 
the nonrelativistic approximation for these observables setting the stage for the complete 
NRQCD analysis. In the last subsection we dwell upon the peculiarities of introducing the 
finite width of the top quark. 

2.1 The effective theory description of the heavy quark threshold 
dynamics. 

Near the threshold the heavy quarks are nonrelativistic so that one may consider both the 
strong coupling constant and heavy quark velocity as small parameters. The threshold 
expansion of the QCD loop integrals has been developed in [|T^. However, to take into 
account the singular threshold effects properly one has to go beyond the finite order QCD 
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perturbation theory. For this purpose the expansion in (3 should be performed directly in the 
QCD Lagrangian within the effective field theory framework. The first step to construct the 
effective theory is to identify all the scales present in the problem. The threshold dynamics 
is characterized by four different scales 0]: 

(i) the hard scale (energy and momentum scale like m^); 

(ii) the soft scale (energy and momentum scale like /5mq); 

(iii) the potential scale (the energy scales like while the momentum scales like (3mq); 

(iv) the ultrasoft scale (both energy and momentum scale like The ultrasoft scale is 
only relevant for gluons. 

By integrating out the hard scale of QCD one arrives at the effective theory of NRQCD 
P). Because the NRQCD Lagrangian does not contain explicitly the heavy quark velocity 
the power counting rules are necessary to construct the regular expansion in this parameter. 
The list of the power counting rules for dimensionally regularized NRQCD and their relation 
to the threshold expansion [jl4| can be found in |]l6l Integrating out the soft modes and 
the potential gluons of NRQCD one obtains the effective theory of potential NRQCD ||13 



which contains potential quarks and ultrasoft gluons as active particles and is relevant for 
the analysis of the threshold effects. In potential NRQCD the dynamics of the quarks is 
governed by the effective nonrelativistic Schrodinger equation and by their interaction with 
the ultrasoft gluons. To obtain a regular perturbative expansion in (3 this interaction should 
be expanded in multipoles. Note that in the process of scale separation some spurious 
infrared and ultraviolet divergences may appear at intermediate steps of calculation which 
cancel each other in the final results for physical observables. The dimensional regularization 
has been recognized as a powerful tool to deal with these divergences [|], |12|, |13|, |14|, |15|, |16 



32|, g3|, ^, g^. 

If the ultrasoft effects are neglected the propagation of a quark- ant iquark pair in the 
color singlet state is described in the potential NRQCD by the Green function (^(x, y, E) of 
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the Schrodinger equation 

(7^-E)G(x,y,E)=5(x-y) (1) 

where Ti is the effective nonrelativistic Hamiltonian. Near the threshold the singular (a^/ 
Coulomb terms should be summed up in all orders in a^. Thus, in threshold region one has 
to develop the expansion in (3 and as around some solution which incorporates properly the 
threshold effects, for example, around the nonrelativistic Coulomb solution. In this case 
the leading order approximation for the nonrelativistic Green function is obtained with the 
Coulomb Hamiltonian 

nc = -— + Vc{x) 

nit 

where Ax = d^. is the kinetic energy operator and Vc{x) = —Cpoislx is the Coulomb po- 
tential, X = |x|. The harder scales contributions are represented by the higher-dimension 
operators in Ti and by the Wilson coefficients of the operators of the nonrelativistic Hamil- 
tonian leading to an expansion in (3 and as- On the other hand the radiation/absorption of 
the ultrasoft gluons by the interacting quark- antiquark pair, the effect of retardation, does 
not contribute in NLO and NNLO (the leading ultrasoft effects in heavy quarkonium have 
been recently computed in ref. [0). Thus, the nonrelativistic Green function of eq. (|1|) is 
the basic object in NRQCD analysis of the threshold effects up to NNLO. In Sections 2.2-2.4 
we relate the observables of e"'"e^ — > ti annihilation and 77 — > ti collisions in the threshold 
region to this Green function. 

2.2 Cross sections. 

We study the normalized cross sections of the top quark-antiquark pair production in e^e~ 
annihilation 



cr(e+e~ —>■ fi+fi~ 
and in 77 collisions 

(t(77 ^ tt) 
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where the lepton cross section 



3s 

is the standard normahzation factor with Uqed being the fine structure constant. Here y/s is 
the total energy of the colhding particles (electrons or photons) in the center of mass frame. 
For unpolarized initial states the following decomposition of the total cross sections is useful 

R-{s) = -^B!'{s) + DaR'{s), (2) 

W(s) = ^""Mt^^'M (3) 

where IT (i?") corresponds to the top quark vector (axial) coupling in e^e~ annihilation 
while R^^ (IV ) corresponds to the colliding photons of the same (opposite) helicity in 77 
collisions. Dy^A are the standard combinations of electroweak coupling constants (see below), 
qt is the top quark electric charge. The cross section for the polarized electron/positron initial 
states is discussed, e.g. in ref. . 

Near the threshold the cross sections are determined by the imaginary part of the corre- 
lators of the nonrelativistic vector/axial quark currents which can be related to the nonrel- 
ativistic Green function and its derivatives at the origin. In NNLO the (potential) NRQCD 
provides one with the following representation of the cross sections 

^ = I c-{a,) + ^^(0, 0, fc), (4) 

mi \ mi) 

47r/V 

= ^C'^ia,)dl^lmGi^, y, k) (5) 

R++^s) = ^l!^^^ ( ( C++ (as) + B++^] lmG{0, 0, k) + dl^lmGi^, y, k)] (6) 
m^ W m^ J J 

R+-^s) = ^^^^C+-ias)dl^lmGi^,y,k) (7) 

where k'^ = —mtE, E = ^/s — 2mt is the energy of a quark pair counted from the threshold 

2mt. A symbolic notation d^y is used for the operator 

3 

i=l 



that singles out the P partial wave of the Green function. The standard electroweak factors 
read 

Dv = QeQt + '^qeqtVeVtd + {vj + al)v^d'^, 
Da = {vl + al)ald^, 

Qe = —1, fe = — 1 + 4 Sin^ 9iv, tte = —1, 

= f , Vt = l- 8 sin^ 9w, at = 1, 



16 sin^ Ow cos2 Ow s - 



The coefficients C*(as) and are the parameters of NRQCD which are responsible for 
matching the effective and full theory cross sections in the limit of weak coupling in NNLO. 
The coefficients 

C\as) = 1 + c\Cf— + e^Cp + . . . 

account for the hard QCD corrections and are determined by the corresponding amplitudes 
with on-shell heavy quarks at rest. The numerical values of these hard coefficients in the 



NLO approximation have been known since long ago [37, 38, 40|. They are explicitly 



= -4, c\ = -2, 



The coefficient has recently been computed in NNLO in different schemes [^2|, |26|, ^| 



Starting from NNLO the hard coefficients acquire the anomalous dimensions and the calcu- 
lation of the NNLO correction requires an accurate separation of hard and soft contributions. 
At the same time these coefficients do not depend on the normalization point of the strong 
coupling constant in NNLO and one can use the different normalization points of as entering 
in the coefficients C* (the "hard" scale fih) and the nonrelativistic Green function (the "soft" 
scale fis), see Section 3.1. 

The coefficients in eqs. (^ and (|^) describe the pure relativistic corrections to the 
cross section which appear when the cross section is evaluated in terms of the correlator 
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of nonrelativistic quark currents. Because the corresponding correction first appears in the 
order 0(/?^) the coefficients can be taken in the leading order in a^. The coefficient is 
related to the nonrelativistic expansion of the vector current and is known P] to be equal to 

4 

B"" = -. 
3 

The calculation of the coefficient B^~^ necessary for the consistent description of the 77 
cross section within NRQCD in NNLO is more involved because the amplitude of 77 ti 
transition is determined by the nonrelativistic expansion of a T product of two vector currents 
?T| , This coefficient, however, can be found by direct comparison with the relativistic 



expression for the cross section expanded in the velocity of the heavy quark (see Section 3.1). 
For the noninteracting quarks (the Born approximation) one obtains the following results 



for the cross sections {(3 = y^l — Amf/s) 

i?++(/5) = QqtN,eiP')i(3 + Oi(3')), i?+"(/3) = 8qtN,eiP')i(3' + 0(/5^)). 

Note that the cross sections K" and R~^~^ are saturated with the S wave contribution and are 
proportional to the Green function at the origin while R"^ and R^~ parts are saturated with 
the P wave contribution and are proportional to the derivative of the Green function at the 
origin. As a consequence they are suppressed in comparison with R" and R'^^ by factor (3"^. 
In the present paper we study the corrections to the total cross sections i?"^ and R^ up to the 
NNLO of NRQCD. Thus i?" is a NNLO contribution to the total cross section i?^ and only 
the leading contribution to R°' is important. On the contrary, the R^~ part can be separated 
from R'' by fixing the opposite helicities of the colliding photons. This makes possible the 
direct study of the P wave production and, therefore, the evaluation of the corrections to 
i?"*"" cross section is of practical interest. 

Concluding this subsection we should also mention that the electroweak corrections to 
the cross sections are known to the one- loop accuracy. For e~^e~ annihilation they have been 
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obtained in ref. and for 77 collisions in ref. ^ . 



2.3 Forward-backward asymmetry. 

The important parameter related to threshold production is a space asymmetry of the dif- 
ferential cross sections. This parameter gives more detailed information on the process and 
allows one to obtain independent experimental data for further testing the theory. The 
forward-backward asymmetry of the top quark production is defined as a difference of cross 
sections averaged over forward and backward semispheres eith respect to the electron beam 
direction devided by the total cross section. A nonvanishing asymmetry appears in e+e~ 
annihilation due to the axial coupling of the top quark to the Z-boson. The expression of 
this parameter for energies near the threshold is given by 

A.B = ^(i + '^^)m (8) 

where 

Eva = qeqtaeatd + 2veaeVtat(f 

is the electroweak factor. The expression for the asymmetry in eq. is given in NLO and 
the explicit correction of order as is taken in the linear approximation that leads to the 
manifest difference of axial and vector hard coefficients in this order. 
The dynamical quantity is a function 

^ 1 ReJG*{p,k)F{p,k)p^dp 

mt JG*{p,k)G{p,k)p^dp ^ ^ 

that describes the overlap of the S and P partial waves. Here pF{p, k) and G{p, k) are the 
Fourier transforms of i9yG(x, y, k)\y=Q and G{x, 0, k) correspondingly. In the Born approxi- 
mation the expression for the function $(/?) can be found in the explicit form and is rather 
simple = Re/3. It vanishes for the real values of energy below threshold. 
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2.4 Top quark polarization. 

The longitudinal top quark polarization in the process e^e~ — > tt averaged over the produc- 



tion angle reads H, H 



where 

P>VA = QeqtVeatd + (f ^ + a^^)vtat(f 

and ^{k) is given by eq. This function enters also the expression for the averaged 

longitudinal top quark polarization in 77 tt process with the same helicity colliding 
photons 

where + (— ) correspond to the positive (negative) helicity photons. 

The extension of the above expressions to the general electron/positron polarization and 
photon helicity and to other component of the polarization vector can be found in the 
literature (e.g. refs. 0, 0). 



2.5 The effects of the finite top quark width. 

As has already been mentioned the sufficiently large t-quark decay width suppresses the 
nonperturbative effects of strong interactions at large (~ 1/Aqc^i) distances and makes the 
perturbation theory applicable for the description of the t-quark threshold dynamics. The 
near-threshold dynamics is nonrelativistic and is rather insensitive to the hard momentum 
details of t-quark decays. Therefore as the leading order approximation the instability of 
the top quark can be parameterized with the constant mass operator. The finite top quark 
width can then be taken into account by the direct replacement nit ^ nit — i^t/'^ in the rel- 
evant argument s — Amf describing the functional dependence of physical quantities near the 
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threshold, ot E ^ E + iTt This approximation accounts for the leading imaginary elec- 
troweak contribution to the leading order NRQCD Lagrangian. Since the essential features 
of the physical situation are caught within this approximation we neglect the electroweak 
effects in higher orders in the strong coupling constant and heavy quark velocity. However, 
in the case of P wave production and S — P wave interference the above prescription is not 
sufficient for a proper description of the entire effect of the non-zero top quark width and 
more thorough analysis is necessary (see Sections 3.2, 3.3). 

In the context of the top quark finite lifetime we should also mention the unfactorizable 
corrections due to the top quark interaction with the decay products which are suppressed 
in the total cross sections |^ but should be taken into account as NLO contributions to the 



angular distribution and top quark polarization [46 



3 Nonrelativistic Green function beyond the leading 
order. 

The basic quantity in the analysis of the threshold effects is the nonrelativistic Green function 
of the Schrodinger equation (|l]). The Green function has a standard partial wave decompo- 
sition 

oo 

^(x, y, k) = Y,{21 + l)ixyyPii^y/xy)Gi{x, y, k) (10) 
1=0 

where Pi{z) is a Legendre polynomial. The partial waves of the Green function of the pure 
Coulomb Schrodinger equation G'^(x, y. A;) are known explicitly 

271 ^0 (m + / + 1 - z/)(m + 2/ + 1)! 

where u = X/k, X = asCpmt/2 with is taken at the soft scale iig- L^{z) is a Laguerre 
polynomial which is chosen in the form 
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We, however, need to know the nonrelativistic Green function for the NNLO Hamiltonian of 
the following form 

n = Hc + AH. 

The second term of this expression describes the corrections to the Coulomb Hamiltonian 

An = + A.Vix) + A2Vix) + A^AVix) + ABpVi^, d^, S). (12) 

The first term of the equation is the standard correction to the kinetic energy operator, 
AnaV{x) = —CACpal/{2mtX^) is the so called non-Abelian potential of quark-antiquark 
interaction and A5j7V^(x, 9x, S) is a standard Breit-Fermi potential known since long 



ago (only the overall color factor Ci? is new). The Breit-Fermi potential contains the quark 



spin operator S, e.g. |^8[. In NNLO the cross section K" is saturated by the final state 
configuration of tt pair with I = 0, S = 1 while R'^^ cross section is saturated by / = 0, 
S = configuration. The Breit-Fermi potential takes the following form when considered on 
the / = states 

AbfV{x) = J + A 5— 5(x) 

X mi mf 

where A'" = 11/3 corresponds to the spin one final state of e~^e^ — > ti production and 
^4++ = 1 corresponds to the spin zero final state of 77 ti production. 

The terms AiV {i = 1, 2) represent the first and second order perturbative QCD correc- 
tions to the Coulomb potential 

A^V{x) = ^Vc{x)iC'o+Cl\n{xf^s)), 

A2V{x) = [^^yVc{x){C',+C',Hx^^,) + C|\n\x^l,)) 



where 



C'o = ai + 2(3oiE, C\ = 2f3o 



Co = [y + J Po + 2(/3i + 2PoaihE + as, 
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C! = 2(/?i + 2/?oai) + S/S^E, Cl = A^l 
31^ 20 

/4343 . TT^ 22 ^2 /1798 56 



V 162 

- (f - 16C(3)) C^TEUf + (yr^n/)' , 
/5i = - y^aTfUj - ACpTFUf. 

Here is defined in MS renormalization scheme. The invariants of the color symmetry 
SU{?>) group have the following numerical values for QCD: Ca = 3, Cp = 4/3, Tp = 1/2, 
77,/ = 5 is the number of light quark flavors, /5o = 11Ca/3 — ATpnf/S is the flrst /3-function 
coefficient, 7^; = 0.577216 ... is the Euler constant and ({z) is the Riemann (^-function. The 
solution to eq. ([^) with the Hamiltonian (|12D can be found within the standard nonrelativistic 



perturbation theory around the Coulomb Green function as a leading order approximation 

^(x, y, k) = Gc(x, y, k) + AG(x, y, k), 

AG'(x, y,k) = -J Gci^, z, k) AH Gc{z, y, k)dz + ... (13) 

In the previous section the threshold effects in the basic observables were reduced to three 
universal functions: the Green function at the origin which is saturated by the S wave 
contribution, the derivative of the Green function at the origin which is saturated by the P 
wave contribution and the function $(/c) which describes the S — P wave interference. These 
functions are analyzed in detail in Sections 3.1-3.3 

3.1 S wave production. 

Only the / = component of the Green function (|1^) contributes to its value at the origin 

G(0,0,A;) = Go(0,0,A;). 
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The explicit expression for the Coulomb part of the Green function has the form 
^^(x, 0, k)U^o = ^ - 2A In {2xfxj) - 2A (^|^ + In 



+ 27e-1 + ^i(1-z/) jj (14) 

where \&n(^) = d"'\nr{z)/dz^ and T{z) is the Euler F-function. The energy independent 
finite part of this expression is chosen for later convenience. Eq. ([l^ ) can be most easily 
obtained from the general expression for the Coulomb partial waves 

C^ix, 0, k) = ^(2A;)2'e-'=^r(/ + 1 - u)U{l + 1 - z/, 2/ + 2, 2kx) (15) 

where U{a, b, z) is the confluent hypergeometric function. In the short distance limit x 
the Coulomb Green function G''^(x, 0, fc) = G^(x, 0,/c) has 1/x and ln(x) divergent terms. 
These terms, however, are energy independent and do not contribute to the cross section. 
Hence these terms can be subtracted without affecting any physical results. The quantity 
fif in eq. (|14D is an auxiliary parameter, the factorization scale, which drops out from the 
physical observables. 

The NLO correction AiG to eq. ([T^) due to the first iteration of AiV^ term of the QCD 
potential has been found in ref. [|r^ where the simple and efficient framework for computation 
of higher orders was formulated. The result of the evaluation of the NLO correction is 



AiGo(0, 0, A;) = 5: F{mf{m + 1) {L,{k) + ^,[m + 2)) - F{m) 



\m=0 m=l n=0 

xF(n)^— + 2 ^ F{m) {L,{k) - 2^e - ^i(m + 1)) - 7sLi(A;) + -U{kf 

where 
and 

Fim) 



m + 1) {m + 1 — v) 
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The NNLO correction A^2^G due to A2V part of the potential and the correction A2'^G due 
to the second iteration of AiV^ part of the correction to the Coulomb static potential have 



been obtained in refs. |]T8|, [T9[. While the technique is rather straightforward the results of 
the calculations are cumbersome and explicit formulae are relegated to Appendix A. 

The method of calculation of the correction to the Green function at the origin due to 
logarithmic terms in the potential is described in details in ref. [^. It is based on the repre- 
sentation of the Coulomb Green function as an expansion over the Laguerre polynomials (p!T|). 
This representation is very close to the standard physical expansion over the eigenfunctions 
that makes the technique transparent and easily interpretable in physical terms. It is equally 
suitable for any partial wave contribution as has been shown in ref. where results for 
the P wave production were found. The results for the S wave part of the corrections were 
reproduced within a different technical framework based on an integral representation of the 
Coulomb Green function in ref. [ ^ . 

The corrections to the Coulomb Green function at the origin due to A^, Vna and Vbf 
terms have been presented in |26, 27]. They are of the following explicit form 



^^Jh+,xL(^^ ^ 



Aa^na,bfG = -j^— ( -A; + 4A ( In ( — I + + ^1 (1 - z/) I - — CFasZ^^2 (1 - i^) 



+ (5-A^ + 2^) Gc{0, 0, . (16) 



In the course of evaluation of this correction to the nonrelativistic Green function one en- 



counters the ultraviolet divergence in the imaginary part of the Green function contained in 
the last term of eq. ([I6|). This divergence is related to the singular behavior of the Coulomb 
Green function at the origin. The particular form of this divergence depends on the reg- 
ularization procedure. The divergence appears in the process of scale separation and is a 
consequence of the fact that the nonrelativistic approximation is not adequate for the de- 
scription of the short distance effects. The hard coefficient C^'^^ computed within the same 
regularization procedure as the Green function itself must have the infrared singular term 
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which exactly cancels the one appearing in the Green function. The hard coefficient can be 
evaluated by matching the effective and full theory cross sections in the weak coupling limit 



17, 26| or by an explicit splitting of the hard and soft contributions using, for example, the 



scale factorization in dimensional regularization [p2| , |3J, ^ . Let us consider the cancelation 
of the divergences and determination of the hard coefficient in the matching scheme. The 
natural regularization for the analysis of the hard part of the corrections is the dimensional 
one |3^, In 4 — 2e dimensions the infrared divergence of the hard contribution in NNLO 



has the form of the ffist order pole in e. The Coulomb Green function at the origin in eq. (|T^ 
can be regularized in the same way. The dimensionally regularized Coulomb Green function 
at the origin can be defined as follows (see Appendix B) 



G^'--(0, 0, k) = --l\^k + 2X\^--+\nl^—j+^E + *i(l - J + (17) 
Note that the Green function regularized in this way differs from that which has been ob- 



tained in ref. I^Sf within another scheme of dimensional regularization. In contrast to eq. (|14D 
in eq. ([T7|) there is no divergence in the Born approximation. The Green function in this 
approximation is a nonrelativistic free propagator and is proportional to k. The ffist order 
pole in e appears only in the ffist order in a^. The O(a^) singular 1/e term in the imaginary 
part of eq. (|T6|) is proportional to Im(G'c'(0, 0, k)) and, therefore, can be absorbed by the re- 



definition of the hard coefficient C^'"^^. For the Green function this redefinition results in the 
substitution G^''' (0, 0, /c) GQ{0,0,k) in eq. (p!6D where the "subtracted" Green function 
reads 

G^(0,0,fc) = G^^'-(0,0,A;)-^l. (18) 

Within the redefined hard coefficient the O(a^) "ultraviolet" 1/e term stemming from the 
corrections to the Green function eq. (pIBf ) exactly cancels the O(a^) "infrared" 1/e term of 
the hard part of the corrections regularized in the same way. Then the (finite) coefficient 
C^'~^^ can be found directly by matching the effective theory expression for the cross sections 
and the result of perturbative QCD calculation of the spectral density in the formal limit 
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cts -C /? ^ 1 up to the order for fih = fJ^s- Eqs. (0, 0) in the matching hmit take the form 
= In^Q^P ((l + (1 - P') + C,^ (y ^ + cl + ^Z?) 



R++ = QN,qt, 



Wtf^ ( (i 



(19) 



where the terms of the relative order O 



. (20) 
are kept. In eq. (pO|) the P wave contribution 



is due to the derivative term in eq. to be discussed. By comparing eq. (|T9|) with the 
NNLO QCD resuh for the cross section K" expanded in the velocity of the heavy quark near 



the threshold |^2|, ^ one finds 



(21) 



where the coefficient has been obtained in refs. ^6|, ^7 



Si 



'39 , 47r2 , 357r2\ ^ 
--C(3) + -ln2-— C, 



151 13,, , 8tt^ , 1797r2\ ^ 
^ + yC(3) + -ln2-— C. 



/44 47r2 11 \ ^ 

The first logarithm in eq. (^Tj) is determined by the renormalization group running of the 
strong coupling constant in the hard momentum regime and is proportional to the first 
coefficient of the /3 function. Thus, both the hard coefficient and the imaginary part of 
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the Green function do not depend on the normahzation point of in the fixed order of 
perturbation theory so one can use different scales for as in these quantities. The second 
logarithm corresponds to the anomalous dimension of the hard coefficient and precisely 
cancels the factorization scale dependence of the Green function due to eq. (|16D making 
the total result independent of the factorization scale. Note that the use of different hard 
and soft normalization scales leads to the incomplete cancelation of the factorization scale 
dependence which, however, is the higher order (0(af)) effect. 

As it has been mentioned above one can bypass the direct matching by the consistent use 
of the same subtraction scheme within dimensional regular izat ion for both hard coefficient 



and Green function. Then matching is automatic In this approach the hard coef- 

ficient is completely determined by the hard renormalization coefficient of the nonrelativistic 
vector current [^. However, in order to compute the corrections to the Green function in 
this case one has to define accurately the Breit-Fermi Hamiltonian and the Green function in 



3 — 2e dimensions P2| , |35(| (in our analysis we use three dimensional Breit-Fermi Hamiltonian 
and direct regularization of the Green function therefore the matching is necessary in order 
to fix the constant relating the two regularization schemes). 

The NNLO analysis of the R^^ cross section is still absent and the constant in the hard 
coefficient is unknown. The logarithmic part of the NNLO contribution to C++(as) reads 

4+ = 4+ - c++^ In f + vr^ {2C, + C^) In l"^] (22) 

where 4^ is a constant to be determined. The relativistic correction to this cross section, 
however, can be extracted from the calculations presented earlier in the literature. Compar- 



ing the known result 



i?++(/?) = QqtN^P + + 0{(3') ) (23) 



with our expressions from eqs. ( pO| , p^ ) we find 

= 1 

3 
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The Green function at the origin can be written in the form which includes only single 
poles in the energy variable. This form seems to be natural for a Green function of the 
nonrelativistic Schrodinger equation 

where ipQm,E'{S^) is the wave function at the origin, the sum goes over the bound states and 
the integral goes over the states of a continuous part of the spectrum. In this way the 
corrections to the Green function stemming from the discrete part of the spectrum reduce 
to corrections to Coulomb bound state energy levels 

and to the values of Coulomb bound state wave functions at the origin 



where 



mt(m+ 1)2' '^O'"' " 7r(m + l)3' 

= ^A.'^,NA,BFEorn + ^2 ''-^Om + ^2 ''-^Om; 

and Afc2^/'Q^ is the correction due to relativistic corrections parameterized by the coefficients 
which we include into the definition of the wave function. 

In NLO an explicit analytical expression for the corrections to the bound state parameters 
has the form |0|, P, Q 



TT 

, 2 (^s(3o 



AiEom = (li{m) + ^i(m + 2) 

TT \ 



^iVom = -Tp^ 3Li(m) + ^i(m + 2) - 2(m + l)^2{m + 1) - 1 - 27^ + 



27r V m + 1. 
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where Li{m) = Li{\/ (m + 1)). The expressions of the NNLO corrections to the energy levels 



2T| , p% and wave functions at the origin [20, 21] are rather cumbersome and given in 
Appendix C, D. 

The continuum contributions in eq. ( pi]) can be directly found by subtracting the discrete 
part of these equations expanded around the Coulomb approximation up to NNLO 



Om 



^{l-E/EUV 

from the result obtained within the nonrelativistic perturbation theory for the Green function 
at the origin (p]3| ) multiplied by (1 — B'^E/rrit). This procedure removes the double and 
triple poles from eq. (|T3|) and leaves only the single poles in the expression for the Green 



function p^). 

An important consequence of the relatively large top quark width is that the most of 
Coulomb resonances are smoothed out. The numerical analysis shows that only the ground 
state resonance in the cross sections is distinguishable. Its separation from others is not 
completely covered by the infrared cutoff provided by the top quark width. Indeed, using 
the pure Coulomb formulas for estimates within the order of magnitude we find 

to be compared with the top quark width Tf = 1.43 GeV. The second spacing between radial 
excitations for the I = partial wave and the first spacing for the / = 1 partial wave are, 
however, much smaller 

li^o^l - Eg\ = \Eg - Eg\ = ^— ^ 0.11 GeV. 

Therefore the contributions of higher resonances are completely smeared out with the top 
quark width. 
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In the limit of vanishing top quark width the NNLO approximation for the energy of the 
first resonance in e~^e~ annihilation reads 

\ 2 / 2 

i^oo = + (Li(A) + 1 - 7i.) + (^) (2C',iL,{\) + 1 - Je) 

+{Cir (^(Li(A) - ^eY + 1 - Y - ^3(1)) + 2C| (^(L(A) + 1 - ^eY - ^ + y)) 

where 



^(A) = ln(|^), L2(A) = ln(^ 



2A ; 

This value is related to the energy of the resonance of the top quark production in 77 collision 
by the hyperfine splitting 

4 A^ 

^00 ~ ^00 — q~'-'F'^s ■ 

6 nit 

The convergence of the perturbation theory series (^) is not fast. For some typical values 
of the soft normalization scale the series for the resonance energy reads 

E^oo = ^0^0^(1 + 0.36 + 0.30), /x, = 25 GeV, 

= E^^{1 + 0.58 + 0.38), /i, = 50 GeV, 
E^o = El^^{l + 0.68 + 0.45), /x, = 75 GeV, (26) 
E^Q = (1 + 0.74 + 0.50), fis = 100 GeV. 
The poor convergence of the series for the resonance energy can be assigned to high infrared 



sensitivity of the pole mass (see, for example, [Q). The convergence can be manifestly 
improved by removing the pole mass from the theoretical expressions in favor of some less 
infrared sensitive mass parameter, for example, the short-distance [EO], potential-subtracted 



or IS" mass |24|. Note that in a finite order of the expansion all the mass definitions are 



perturbatively equivalent. The infrared safe mass parameters, however, are "closer" to the 
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physical observables since in contrast to the pole mass the corresponding perturbative series 
are supposed to be better convergent (less divergent). 

Due to the finite top quark width the location of the peak (maximum) of the cross section 
is not given only by the position of the ground state resonance but is also affected by the 
contribution of the higher (smeared out) resonances and the continuum contribution. Due 
to this effect the absolute value of the NNLO peak energy (^Sj) counted from the threshold 
is less than the absolute value of the energy of the ground state resonance -Eqo^^ by about 
200 MeV i.e. ~ 7%. This shift is essentially smaller than the one related to the perturbative 
QCD corrections to Coulomb values but considerably larger than the leading nonperturba- 
tive contribution due to the gluon condensate ||55| which is suppressed parametrically as 

3.2 P wave production. 

The derivative of the Green function at the origin is saturated with its Z = 1 component and 
explicitly given by the relation 



For the Coulomb Green function from eq. ([T5|) we obtain the closed formula for the partial 



wave / = 1 Green function at the small space separation of particle 

+A (2(.^ - A^) (± + „> (A) + 20. - H + (1 - .)) 4)) (27) 

where /i/ is the analog of the parameter fif for the / = partial wave. In the short distance 
limit X ^ the derivative of the Coulomb Green function (or the partial wave with / = 1) has 
(n = 1,2,3) and ln(x) singularities. In contrast to the case of the S wave production, 
the value at the origin for the P wave partial Green function contains divergent terms 
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that explicitly depends on the energy (or wave vector) k. However, these terms do not 
contribute to the cross section for the vanishing top quark width = because they have 
no discontinuity across the physical cut in the complex plane of the energy variable in the 
approximation of top quark zero width. The case of the non-zero top quark width requires 
to perform a more detailed analysis given below. 

The correction to the / = 1 partial wave at the origin due to the first iteration of AiV 



term of the QCD potential has been found in ref. |^ 



AiGi(0, 0, k) = ( E F{m)\m + l)(m + 2)(m + 3) {L,{k) + ^i{m + 4)) 



\m=0 



-2 f; F{m)F{n) + ^) + ^) + ^) + 2 f] F(m) ( 2Jo(m) + (m + l)(m + 2)Li(A;) 



m=l n=0 ^ ^ m=0 



where 



+ (1 + z/)(Ji(m) + (m + l)L^{k)) + ''-^^^±^{J^{m) + 2L^{k)) ) + I{k) 



[m) 



(m + 2 - z/)(m + l)(m + 2)(m + 3) 
Mm) = -2^i(m + 1) - A^E + 3, 

Ji(m) = (m + l)(-^i(m + 2) - 2-fE + 2), 
~ , , (m + l)(m + 2) ( ^ , , 3 



2 V ^ ' 2 

Kk) = ^ 12 ~ ^ ^^^"^ 4 2 ^^'^ 

9 -127s 2 3 3 1 A r , ^ 2 



+ (^7i. - 1 - 3Z/ + + -u' + -z/^j Li(fc) + (^-- + -z.^ j L^iky 

For the derivative of the Green function at the origin (or for the / = 1 partial wave) the 
analog of eq. (|2^) reads 

where, symbolically, 
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Here Eim is the / = 1 bound state energy. In NLO approximation these quantities read |^ 

^- ^1 + ^ 



\ 2 

^1- = ? ^ + r^2C^ (L^{m + 1) + ^i(m + 4 



and 



I^L(0)r = , ^ 1 + T^C^ 5Li(m + 1) + 5^i(m + 4) - — (m + 2) - 1 



TT (m + 2)5 \^ 47r ' V 3 



m— 1 



_(n + l)(n + 2)(n + 3) 

n=0 



+2 V 

(m + l)(m + 3)(m — n)^^ 



The continuum contribution is obtained in the same way as it was done in the previous 
section for the S wave production. 

In the case of P wave production the simple shift E E ^ iV \\y the nonrelativistic 
approximation is not sufficient to describe properly the entire effect of the non-zero top 
quark width p[. Indeed, eq. p7| ) in the limit x — with the nonvanishing width has 
the divergent imaginary part with the leading power singularity ~ F^/x related to the free 
Green function singularity and the logarithmic singularity ~ F^a,, ln(x) produced by the 
one Coulomb gluon exchange. The presence of these singularities clearly indicates that the 
coefficient of the constant term linear in F^ gets a contribution from the large momentum 
region and cannot be obtained within the pure nonrelativistic approximation. Like the hard 
coefficients it should be computed in relativistic theory. Parametrically this contribution 
is not suppressed in comparison to the pure nonrelativistic contribution in the threshold 
region. At = 0, for example, the ratio between the relativistic (proportional to F^) and 
nonrelativistic (Coulomb) contributions is of order ~ 1. Since we are interested in the 
NLO corrections the O(Ftas) term also has to be taken into account. By construction, the 
nonrelativistic effective theory has to reproduce the perturbation theory predictions in the 
formal matching limit a^, /5 -C Fj/m^ ^ 1 where both effective theory and perturbation 
theory descriptions are valid. Thus one has to compute 0(Ft) and ©(F^cts) terms within the 
relativistic perturbation theory and then to fix the parameters of the effective nonrelativistic 
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theory so that it reproduces the perturbative results in the matching hmit. Within the 
relativistic perturbation theory the relevant contributions can be obtained by inserting the 
complex momentum-dependent mass operator into the top quark propagator at /3 = (only 
the leading terms in Tt/rrit should be kept). In the leading order in this procedure has 
been done in refs. [|]. The result reads 

777 ^ 

Gf(0,0,fc)|r, = -^^7ir, 

ODTT 

where gi is a coefficient coming from the relativistic treatment with numerical value gi = 
0.185 . . . For the term of the order 0{rtas) the necessary calculation has been performed in 
ref. [^. It has been shown that the proper relativistic analysis leads to fixing the auxiliary 
parameter of eq. p7| ) /i/ = (72^* where g2 is the coefficient coming from the relativistic 



treatment. Its numerical value is g2 = 0.13 . . .. 

Here we should note also the problem of the previous numerical analysis of the P wave 
contribution [Q. While solving the Schrodinger equation ([l|) numerically for the finite 



top quark width one has to introduce an explicit ultraviolet cutoff for the nonrelativistic 
expressions divergent in the large momentum region. To get rid of the cutoff dependence 
one has to compute the hard contribution within the relativistic approximation using the 
similar prescription for the infrared cutoff. This, however, has not been done and as a 
consequence the 0{rt) and 0{rtas) contributions to the cross section were not determined 
within the numerical framework of ref. On the other hand the total 0{rt) contribution 



to the cross section is numerically small in comparison with that of the regular completely 
nonrelativistic terms of eq. (|27| ) which saturate the total result for the energies not far below 
the threshold. 



^Recently the 0{Tt) contribution has been estimated within the numerical approach ||30| by using the 
physical (relativistic) phase space for the unstable top quark to regularize the divergence of the nonrelativistic 
approximation. 
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3.3 S — P interference. 

In the zero width approximation the function (|^) allows for the following decomposition 

<I>{E) = %,i{E) + 

where $con and ^poi are the continuum and bound state poles contributions correspondingly. 



This is known [0, Q that the continuum contribution is not affected by the Coulomb effects 



and above the threshold one has the Born approximation result 

^Z^iE) = Re./^ 

V m 

even for the Coulomb Green function in eq. (Q). Below the threshold in the Coulomb ap- 
proximation one gets 

where the quantities 

c _ (m + l)(m + 3) 
'^"^ ~ mtTi (m + 2)5 

measure the overlap of the S and P wave functions. Note that in the zero width limit the 
function ^^^i does not vanish due to the Coulomb degeneration of the energy levels with 
different /: -E^+i = -E'lm- It was pointed in ref. that the continuum contribution gets 
no soft corrections. Thus in NLO for a finite top quark width we have the simple result 



<^>con{E) = Rex 



lE + zTt 



mt 

The corrections to the pole contribution are less trivial. They can be computed using the 



powerful technique developed in refs. ||T^, |21], |31| . The result reads 



CP,„H^J Ke |^2L. ^E,^^^^^ -E + ^T,){E,^ -E- ^T,) ) [i^, {E,^ - Ef + Tl) ^"'^ 
where 



1 + 4L(m + 1) + 4^(m + 4) + 2 (m + 2) 

vr V ^ ' ^ ' m + 3 6^ ' 
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degeneration is lifted by the logarithmic corrections to the potential. Strictly speaking, our 
approach is valid only if the level splitting Eom+i — -Eim is much smaller than the top quark 
width (which is realized for the actual values of these quantities). For Eq^+i — -Eim > 
the nonrelativistic analysis is not applicable for the S — P interference below the threshold 
because the double poles of eq. (|28D disappear and the nonrelativistic contribution is not 
enhanced in comparison to the relativistic one in this case. 

Note that for the finite top quark width the interference of the free I = and / = 1 
partial waves results in the logarithmically divergent 0(Tt) term in the numerator of eq. (P) 
(this term does not include the factor Tt explicitly but it is suppressed in comparison to 
the leading term which is proportional to l/F^ as the leading term in the denominator of 
eq. (1^) for the free quark Green function). This divergent term is of the same nature as 
the divergence in the P wave amplitude discussed in Section 3.2. An accurate calculation 
of this term can be done only within the relativistic approximation. In contrast to the P 
wave production this term is parametrically suppressed above the threshold in comparison 



to the nonrelativistic continuum contribution at least by the factor yFt/m^ at E ~ and 
can be safely omitted. However, it becomes important below the resonance region when the 
nonrelativistic contribution becomes small. Moreover, the denominator in the right hand side 
of eq. (^ decreases rapidly below the ground state pole. Therefore a small uncertainty in 
the numerator would lead to a large uncertainty in the function $(A;) and a reliable estimate 
of its numerical value is not possible in this region within the nonrelativistic approximation. 
Strictly speaking the accurate determination of the function $ below the ground state pole 
requires the calculation of the relativistic 0(Tt) contribution to the S wave cross section 
(the denominator of eq. (^) which is not usually considered since it does not lead to the 
divergence in the nonrelativistic expression. 
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4 Discussion. 



The results of the numerical analysis for the physical observables based on the obtained 
analytical expressions are plotted in Figs. 1-4. 

The constant c^^ appearing in the hard coefficient C^^{as) in the O(a^) order remains 
unknown. The calculation of this parameter is necessary for the formal completion of the 
NNLO analysis. To find its numerical value one has to compute the 0{a1) perturbative 
QCD correction to the 77 cross section near the threshold in the formal limit ^ /5 ^ 1 
and compare it with 0{a'j.) term in eq. (|20|) . In the case of e~^e~ annihilation, however, 
the analogous contribution parameterized by is relatively small (about 10% of the total 
NNLO correction) and the correction to the physical observables in NNLO is saturated with 
the soft part of the total contribution determined by the corrections to the parameters of 
the nonrelativistic Green function. Thus one can reasonably hope that the similar situation 
can also take place for 77 collisions. However the importance of this parameter for physical 
observables is not crucial, it affects only the overall normalization of the cross sections. For 
example, it does not shift the position of the resonance which is an important characteristic 
of the production and does not enter the ratio R'^'^ (E) / R'^'^ (0) . For the numerical analysis 
of the cross section R'^ we set = 0. 

In our approach we deal with the soft corrections by summing them into the energy 
denominators of the discrete part of the Green function. In other words we treat the soft 
corrections as effective corrections to the parameters of the Green function written in a fixed 
functional form. The same approach has been advocated in refs. p2| , [27| , p9| , |30| , [56[| where all 
the corrections to the Green function have been found (numerically or analytically) in the 



form (|2^) . In ref. however, a part of the NNLO corrections has not been resummed to 
the energy denominators of the discrete part of the Green function. On the other hand, in 
refs. p7| , pOj |56[| the Schrodinger equation (|l]) has been solved numerically, i.e. the NLO 
and NNLO correction to the Coulomb Hamiltonian have been taken into account effectively 
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in all orders of the nonrelativistic series (0) for the Green function while we work strictly 
in NNLO. Our formulae reproduce the numerical result for K" of the most recent numerical 



analysis [|^, ^ ^ with 1% — 3% accuracy that can be assigned to the contribution of 
the higher iterations of the NLO and NNLO corrections to the potential in eq. ([13D beyond 
NNLO. 

For the total cross sections which are dominated by S wave contribution we find the 
typical size of the NNLO corrections to be of the order of 20% in the overall normalization 
of the cross sections and ~ 40% in the resonance energies expressed in terms of the top 
quark pole mass, i.e. of the order of the NLO ones (see Fig. 1). Though the inclusion of 
the NLO corrections leads to a considerable stabilization of the theoretical results for the 
cross sections against changing the normalization point, the NNLO corrections do not lead 
to better stability as compared to NLO. In the overall normalization of the cross sections 
the NLO and NNLO corrections cancel each other to a large extent while the NLO and 
NNLO corrections to the resonance energies are of the same sign and shift the resonance 
farther from the threshold. They also make the peak more distinguishable which is the main 
difference between the leading Coulomb and NNLO approximations. 

The leading order approximation for and W cross sections are the same up to the 
normalization factor 2ql. Up to the overall factor the difference between the cross sections 
is determined by NNLO QCD and relativistic corrections (see Fig. 2). Above the threshold 
this difference is determined by the difference between 5++ and coefficients and between 
P wave contributions to eqs. (0, 0) i.e. by the pure relativistic corrections. Below the 
threshold in the resonance region this difference is determined also by coefficients and is 
quite sensitive to the value of Ug. 

Though using an infrared safe mass parameter instead of the pole mass improves the 
convergence of the series for the resonance energies it does not affect the huge NNLO correc- 
tions to the cross sections normalization. Moreover, it is not clear if there exist physically 
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motivated mass and strong coupling parameters providing fast uniform convergence of the 
perturbative expansion for the cross sections in the threshold region. The absence of such 
a parameterization would mean the unavoidable significance of the high orders terms of the 
threshold expansion. Some high order effects have been already considered in the literature. 
The leading logarithmic corrections related to the renormalization group evolution of the 
hard coefficient C" have been computed ||2^. The corresponding corrections to the R" cross 
section are ±5%. In ref. the running of the strong coupling constant has been taken 
into account by introducing the energy dependent soft normalization point of entering the 
Coulomb potential in the numerical solution of the Schrodinger equation. The resummation 
of the renormalization group logarithms has an essential (up to 10%) effect in the resonance 
region and reduces the normalization scale dependence of the result. Furthermore, the ef- 
fect of retardation which introduces a new type of contributions absent in NLO and NNLO 
has been analyzed for the low lying resonances [^. The characteristic scale of the leading 
ultrasoft contribution was found to be about —5% for the square of the ground state wave 
function at the origin and +100 MeV for the ground state pole position. 

The result for the axial coupling contribution to e"'"e" tt cross section is in a good 
agreement with the numerical analysis of ref. [Q. Up to the trivial normalization this 
contribution coincides with the cross section (Fig. 3). Numerically it does not exceed 
2% of the total cross section and less than the uncertainty due to the normalization scale 
dependence. 

The cross section and the function $(A;) obtain no contribution from the ground 

state resonance and, therefore they are rather smooth because the top quark width smears 
the higher resonance contributions very efficiently (Figs. 3, 4). These quantities is rather 
insensitive to a variation of the normalization scale. A typical NLO correction to 
is about 10% while the one to $(A;) is about 15% (the corrections to the forward-backward 
asymmetry and top quark polarization include also the hard normalization factors which have 
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not been included to itself and the nonfactorizable corrections corrections discussed in 
Section 2.4). Our result for the function ^{k) (Fig. 4) is in a good agreement with the results 
of numerical analysis P, § for the energies above the ground state resonance. There is some 
discrepancy between the results below the resonance. However the reliable estimate of the 
function $ is not possible in this region with the pure nonrelativistic treatment of the top 
quark width as has been explained in Section 3.3. 

The final remark of this section concerns the optimal choice of the normalization and 
factorization scales. The hard scale appears in the hard coefficients as In {nit/ fi^) ^-c the 
typical hard scale of the problem is the top quark mass. Though in a fixed order of the 
perturbative expansion the hard coefficients do not depend on fih one can put fih ~ mt 
to minimize the potentially large logarithmic contributions of the higher order terms. In 
practice the NNLO results are almost independent of fih when fih ~ "^t- On the other hand 
the requirement of convergence of the perturbative expansion around the Coulomb Green 
function restricts the allowed range for the choice of a soft normalization point which can 
be used for reliable estimates. The soft physical scale of the problem is determined by the 
natural infrared cutoff related to the top quark width y/rritTt that measures the distance 
to the nearest singularity in the complex energy plane and/or by the characteristic scale 
of the Coulomb problem A i.e. fis ~ 15 GeV. Both scales are rather close to each other 
for the case of top quark that makes possible a uniform description of both perturbative 
QCD and Coulomb resonance effects. Indeed, for fig ~ 15 GeV the soft NLO correction, 
for example, to the energy level (^) reaches its minimal magnitude. However, at this scale 
the NNLO correction exceeds the NLO one and the series for the energy levels seems to 
diverge. Moreover, for such a low soft normalization point the NNLO corrections to the 
wave function at the origin which cannot be dumped by the quark mass redefinition become 
uncontrollable. This is not surprising since the normalization scale is defined in a rather 
artificial MS scheme that has httle to do with peculiarities of tt physics and there is no 
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reason for a literal coincidence of j^s parameter with any physical scale of the process. The 
relative weight of the NNLO correction to the Green function as well as the dependence of 
the cross sections on fig is stabilized at /is^40 GeV which can be considered as an optimal 
choice of the soft normalization point. The price one pays for using different soft and hard 
normalization scales is the incomplete cancelation of the factorization scale dependence but 
this effect is suppressed by an additional power of a^. Another source of the dependence on 
the factorization scale is the factorized form P) of the cross sections where some higher 
order yiiy-dependent terms are kept. The numerical analysis, however, shows that the results 
are rather insensitive to the factorization scale chosen in the region /ij ~ mt. 

5 Conclusion 

The basic observables of the top quark pair production in e~^e~ annihilation and 77 colli- 
sion have been considered in the threshold region. The threshold effects are described by 
three universal functions related to the S, P wave production and S — P wave interference 
which have been computed analytically within (potential) NRQCD. An explicit analytical 
expression for the soft part of the NNLO corrections to the total cross section has been 
obtained. The e~^e~ tt threshold cross section has been obtained in NNLO in the closed 
form including the contribution due to the top quark axial coupling. The forward-backward 
asymmetry of the quark-antiquark pair production in e^e~ annihilation and top quark po- 
larization in both processes have been computed analytically up to NLO. The running of the 
strong coupling constant and the finite top quark width effects in the P wave production 
and S — P wave interference have been taken into account properly within the analytical 
approach. 

In combination, these uncorrelated observables form an efficient tool for investigating 
quark interactions. As independent sources they can also be used for determination of the 
theoretical uncertainty in the numerical values of the strong coupling constant a<j, the top 
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quark mass, and the top quark width extracted from the experimental date on top-antitop 
production. 

The high order corrections turn out to be relatively large for all observables and important 
for the accurate description of the top quark physics near production threshold. 
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Appendix. 

A. The correction A^2 ^G due to the A2V part of the potential has the form ||18 
A^G = (^)' [t/irnr {{m + 1) {C^, + md + L{kfC\ 

+ (m + l)^i(m + 2) (Cl + 2L{k)Cl) + K{m)Cl 

00 m—l , r) -I- 1 

+2 y y F{m)F{n) -— ( Cl + 1L{k)Ct\ + Kim, n)C 

00 

+2 F{m) (pi + L{k)Cl + [Likf + J{m))Cl - (27^ + ^i(m + 1)) [cl + 2L{k)C, 



m=0 
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+mc'^ + ( -^Em + h{kf ) ci + i{k)c^ 



where 



TT 



2 



ir(m) = (m + 1) ^i(m + 2)^ - ^2("^ + 2) H 



^ 3 (m + l)2; 

-2(^i(m+ 1) +7s), 

+ 1 / . X 1 \ 

Kim, n) = 2 \&i(m - n) h 27^; 

m — n \ n + 1 } 

TTl ~t~ 1 

+2 (^i(m - n + 1) - ^i(m + 1)), 

m — n 

Jim) = 2(^i(m + 1) + 7s)2 + ^2(m + 1) - ^i(m + 1)^ + 27I, 

I{k) = I^^E + y ) m - jEHkf + h{ky. 
The correction A'2^G due to the second iteration of Ail^ term |jl9| 

W) 47r 2/. \to ^ ' 
(Co' + (^(m + 2) + L(A;))C|)' 

/ 00 m— 1 i — 1 ;n _L 1 



Vm=2 /=1 n=0 



(/ — n)(m — n) 



where 



00 m— 1 71—1 ; I 1 



i7(m) = 

m + 1 — 
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B. We define dimensionally regularized value of tlie Coulomb Green function at the origin 
directly through the relation 



with d = "i — 2e. Using the following representation of the momentum space Green function 



8n^Jo V t J + A;2(l + 2t)2)2 

one obtains 



^ ^ ' ' ^ 2-K \ k J \ t J {1 + 2tY^ 
where we omit inessential factors related to the precise definition of integration measure in d 
dimensions. These factors lead to the multiplication of the Green function with an additional 
quantity l + 0{£) and can be taken into account by the redefinition of /xj scale. The integral 
in the right hand side of this equation is 

where B{z,w) is the Euler i?-function and 2Fi{a,b] c; z) is the hypergeometric function. 
Upon expanding the above expression in e around 6 = one arrives at the final result for the 
dimensionally regularized Coulomb Green function. The factorization scale ;U/ in eq. ([T7|) is 
chosen in such a way that it is true as written. Note that the Green function regularized in 
this way does not automatically match the hard coefficient computed in MS scheme of the 



orthodox dimensional regularization |3^ . 



C The NNLO corrections to the square of the Coulomb ^5*1 and ^5*0 heavy quark bound 
state wave function at the origin have the form pll 
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+7e + + 1) - ^ 



{m + l)J J' 

^?^^L = {^C', + L{m)Cl + L{mfCl) + (-1 - 2^^ + ^ + + 2) 



-2(m + l)^2(m + l)){Cl + 2L(m)C2') + ( + 2i^(m) - 2^i(m + 2) 

A?VL = (^)' (3(C„' + (L(m) + *:(m + 2))Ci)^ 

(r m— 1 »^ _|_ 1 

^ m + 1 ^, ,\ ^,^1.2 I ^ n + l ^ m + 1 



^^ig^ (/ + l)(n + l) (/ + l)(m + l) (/ + l)(m + l) 



^ 1^ \(^ — m)^(/ — m)^ (n — m)^(/ — m){n — I) {I — mY{n — m){n — I) ^ 
~ ( (/ + l)(m+l) (/ + l)(m + l)2 



n=m+l 1=0 



[n — my(l — mY {n — mY{l — m){n — l){n + 1) 



(Z + l)(m + l) \+f^ ( (^ + 1)' 



{I - m)^(n -m)(n-l)J ;Sz=tl^i \(^ - - H 

(/ + l)(m + l)2 (m + l)2 



2 



(n — my(l — m)(n — l){n+ 1) (/ — mY{n — m){n — I) ^ 
where L(m) = L{\/{m + 1)) and 

?7(m,n) = 3 + --2- ^ ^ 



m + n + 2 {n — m){n + m + 2) 
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D. The NNLO corrections to the Coulomb ^Si and ^5*0 heavy quark bound state energy- 



levels [20, 21, 53 



^A^,NA,BFEom 



Clal (Ca . 11 1 



(m + 1) \ 



+ 



16 (m+1 



A^E^^ = 2 (Cl + L{m)Cl + L{mfCl + ^i(m + 2){Cl + 2L(m)a 



(m+ 1) 



2 ' 



A5')Eom= (^) ((cl + {L{m)-2 + ^>^{m + 2))Cl) (cl + {L{m) + m ^{m + 2))Cl 



+ 



(m+ 1) 



(7b + ^i(m + 2)) - 2$2(m + 1) - (m + l)^3(m + 1) {Cl 



1\2 



Figure captions 

Fig. 1. The normalized cross section R"{E) in the leading order (solid lines), NLO (bold 
dotted lines) and NNLO (bold sohd lines) for = 175 GeV, Vt = 1.43 GeV, a,(Mz) = 0.118 
and fig = 50 GeV, 75 GeV and 100 GeV. The dotted hne corresponds to the result in Born 
approximation. 

Fig. 2. The normalized cross sections R'^{E) (dotted lines) and W{E) (solid lines) in 
NNLO for ct+ = 0, = 175 GeV, Tt = 1.43 GeV, a,(M^) = 0.118, sin^^^^ = 0.232, 
Mz = 91.2 GeV and /i, = 50 GeV, 75 GeV and 100 GeV. 

Fig. 3. The normalized cross section R~^~{E) in the leading order (dotted lines) and 
NLO (bold solid lines) for rut = 175 GeV, Tt = 1.43 GeV, as{Mz) = 0.118 and /x, = 
50 GeV, 75 GeV and 100 GeV. The solid line corresponds to the result in Born approxima- 
tion. 

Fig. 4. The function $(i?) in the leading order (dotted lines) and NLO (bold solid lines) for 
mt = 175 GeV, = 1.43 GeV, ^.(M^) = 0.118 and /i, = 50 GeV, 75 GeV and 100 GeV. 
The solid line corresponds to the result in Born approximation. 
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